Pairing symmetry in a two-orbital Hubbard model on a square lattice 
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We investigate superconductivity in a two-orbital Hubbard model on a square lattice by applying 
fluctuation exchange approximation. In the present model, the symmetry of the two orbitals are 
assumed to be that of an s orbital. Then, we find that an s-wave spin-triplet orbital-antisymmetric 
state and a p-wave spin-singlet orbital-antisymmetric state appear when Hund's rule coupling is 
large. These states are prohibited in a single-orbital model within states with even frequency 
dependence, but allowed for multi-orbital systems. We also discuss pairing symmetry in other 
models which are equivalent to the two-orbital Hubbard model except for symmetry of orbitals. 
Finally, we show that pairing states with a finite total momentum, even without a magnetic field, 
are possible in a system with two Fermi-surfaces. 

PACS numbers: 74.20.Rp, 74.20.Mn, 71.10.Fd 



I. INTRODUCTION 

It has been recognized that orbital degree of freedom 
plays important roles in determination of physical prop- 
erties, such as colossal magneto-resistance and complex 
ordered phases of manganitesji*^ and exotic magnetism in 
/-electron systemSfSii From experimental and theoretical 
studies of such systems, it is revealed that magnetism in 
multi-orbital systems has a rich variety. In recent years, 
effects of orbital degree of freedom on superconductiv- 
ity have also been discussed theoretically for several raa- 
terialsj^'^ii'^'^'i^iii and it has been found that orbital 
degree of freedom is important, e.g., for determination 
of pairing symmetry in a system. In particular, orbital 
degree of freedom probably plays a role in triplet super- 
conductivity of Sr2Ru04 ,— " in which Fermi surfaces are 
composed of t2g orbitals. 

To understand the effects of orbital degree of freedom 
on superconductivity, it is still important to gain a defi- 
nite knowledge on those in a relatively simple model such 
as a Hubbard model for two orbitals with the same dis- 
persion, although, at present, it is difficult to find di- 
rect relevance of such a simple model to actual materi- 
als. For such a purpose, superconductivity in two-orbital 
Hubbard models has been studied by a mean-field the- 
oryii and by a dynamical mean-field theoryii^ii^ These 
studies have revealed that an s-wave spin-triplet state, 
which satisfies the Pauli principle by composing an or- 
bital state of a pair antisymmetrically, is a candidate for 
a ground state. This fact is in sharp contrast to a single- 
orbital model in which pairing states with even parity in 
the wavenumber space, such as an s-wave state, should 
be spin-singlet ones due to the Pauli principle within 
states with even frequency dependence. Note that odd- 
frequency states are hard to be realized in ordinary cases. 
In addition, in a two-orbital system, it is also possible to 
realize odd-parity spin-singlet states. Thus, the variety 
of pairing states in a two-orbital model is much large. 

However, in the above studies, possibility of supercon- 
ductivity other than the s-wave state is not considered. 
In particular, within the standard dynamical mean-field 



theory, i.e., in infinite spatial dimensions, we cannot deal 
with spatial dependence of a pairing state. Thus it is 
desirable to study superconductivity in a multi-orbital 
model on a finite-dimensional lattice and determine the 
most plausible candidates for pairing symmetry of super- 
conductivity. 

In this paper, we investigate possible superconducting 
states of a two-orbital Hubbard model on a square lat- 
tice by applying fluctuation exchange (FLEX) approxi- 
mation. The FLEX approximation has been extended 
to multi-orbital models i^'^'^ i^^'^"'^ We classify supercon- 
ducting states by spin states, orbital states, and repre- 
sentations of tetragonal symmetry. We also discuss pair- 
ing symmetry in other models which are equivalent to 
the two-orbital Hubbard model. In particular, we find 
that pairing states with a finite total momentum like 
the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) si&te^^^ 
even without a magnetic field, are possible in a system 
with two Fermi-surfaces. 

The organization of this paper is as follows. In Sec. HIl 
we introduce the two-orbital Hubbard model. In Sec. lIIIi 
we explain the FLEX approximation and categorize pair- 
ing symmetry of the model. In Sec. IIVI we show re- 
sults obtained with the FLEX approximation. In Sec. |Vl 
we discuss pairing symmetry in other models which are 
equivalent to the two-orbital Hubbard model except for 
orbital symmetry. We summarize the paper in Sec. IVII 



II. HAMILTONIAN 

To investigate superconductivity in a multi-orbital sys- 
tem, we consider a two-orbital Hubbard model given by 

k.r.cr i,T 

+ U'^ nana + c\^A2<,'(^^^'y'C^2a (j^) 
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where Cira is the annihilation operator of the electron at 
site i with orbital r (= 1 or 2) and spin tr (=t or j), 



is the Fourier transform of it, n,. 



and 



Uir = J2(7 '^ira- The coupliug coustauts U , U' , J, and J' 
denote the intra-orbital Coulomb, inter-orbital Coulomb, 
exchange, and pair-hopping interactions, respectively. In 
the foUowings, we use the relation U = U' + J + J' , which 
is satisfied in several orbital-degenerate models such as a 
model for p-orbitals, a model for Cg orbitals, and a model 
for t2g orbitals We also use the relation J = J' , which 
holds if we can choose wavefunctions of orbitals realfi^ 

Concerning the kinetic energy ekr, we consider only a 
nearest-neighbor hopping integral t for both orbitals for 
simplicity, and the kinetic energy is given by Cki = ek2 = 
Ek — 2i (cos kx + cos ky). Here we have set the lattice con- 
stant unity. We note that we assume that the symmetry 
of orbitals is an s-orbital one, i.e., an orbital state does 
not change by a symmetry operation of a lattice, such as 
inversion. This assumption is crucial to determine pair- 
ing symmetry of superconductivity. 

Here we note that the model Hamiltonian ^ can also 
describe a system with different orbital-symmetry with 
a special condition. For example, a model for doubly 
degenerate Px and Py orbitals with the Slater-Koster in- 
tegrals (ppa) = ippn) = t is given by Eq. ^ with 
Cki = ek2 = Ck- Another model for s orbitals described 
by Eq. ([T]) with eki = — ek2 = Ck is equivalent to the 
model with eki = ek2 = fk, if we change the phases of 
the wavefunctions for t = 2 orbitals by exp[iQ • r^] at 
each site i, where Q — (tt, tt). We also discuss pairing 
symmetry in such equivalent models in Sec. |Vl 



III. FORMULATION 

In this section, we derive equations for response func- 
tions, classify symmetry of superconductivity, and derive 
a gap equation for the anomalous self-energy for each 
symmetry. 

A. Green's function 

First, we derive equations for the Green's function in 
the normal phase. In general, the Green's function is 
defined by 

Gri(Ti;r2<T2 (k, t) = ^(TrCkricri (''')cj^.r2tT2)' 

and the anomalous Green's functions are defined by 

-fricri:r2(T2 (k, t) = — (T7-Ckri(Ti (''")C— kT2(T2 ) i (3) 
PU;r.<r,i^^r) = -{Trcl^^^,^iT)cl^J, (4) 

where T^- denotes the time-ordered product and (• • • ) 
denotes the thermal average. The Heisenberg represen- 
tation for an operator O is defined by 



where iVtot = J2i r is the total number operator of 
electrons and /x is the chemical potential. It is conve- 
nient to use the Fourier transformation with respect to 
imaginary time given by 



0{ier. 



(6) 



where f3 — 1/T with a temperature T and e„ = (2n -|- 
l)7rr is the Matsubara frequency for fermions with an in- 
teger n. Here we have set the Boltzmann constant unity. 
Then, the Dyson-Gorkov equations are given by 

Gra;r'a'{k) ^6rr' 5aa' G'^°\k) 

+G^'Hk)^r.;r"a"{k)Fl„^„._^,^,ik)l 

= J2 [GW(fc)S,,;,„,„(fc)F,,v";r'.'(fc) (8) 
-G(0) {k)cj)^a,r"." {k)Gr'a';r"a" {~k)] , 

where Sro-iT'o-' (k) is the self-energy and ^to-itV (k) is the 
anomalous self-energy. Here we have used the abbrevia- 
tion k — (k, ie„). The non-interacting Green's function 
is given by 



(9) 



In the normal state, the Green's function and self- 
energy do not depend on spin and orbital states, 

i.e., GTa;T'cr'{k) = 6^^,Saa'G{k) and l^rair'a'ik) = 

dTT'Sa-a-''^{k), and then the Dyson-Gorkov equation is 
given by 

G(fc) ^ G^"\k) + G(")(fc)E(fc)G(fc). (10) 
The self-energy is given by 



m^^T.^{q)G{k-q), 



N 



where 



(11) 



(12) 



in the FLEX approximation. Here, N is the number of 
lattice sites and q = (q, iuJm), where ujm = 2nnrT is the 
Matsubara frequency for bosons with an integer m. The 
matrix y normal ^^-j given by 



0(r) = e"(^-''^-')Oe-"(^-''^-'), 



(5) 



^normal = ^[U^'x^iqW ~ W X^'^HqW / 2 + W 



(13) 
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The matrix elements of and If^ are given by 





— ^22:22 




— t^22;22 




(14) 


^11)22 


= ^22;11 








(15) 


"-^ll;22 


— "-'22;11 


= 2U' - 


J, 




(16) 


"-^12;12 


— "^21;21 


= c/', 






(17) 


"-^12;12 


— "-'21;21 




-2J, 




(18) 


TP 

•-^12:21 


— "^21:12 


— '^12;21 


— "-^21:12 


= J', 


(19) 



and zero for the other elements of these matrices. The 
matrices for susceptibilities x^W) for the spin part and 
X'^io) for the charge part are given by 

r{q) = X^°\q)[l-U^X^'>\q)]-\ (20) 

x'{q) = X^'\q)[l + U^X^''\q)]-\ (21) 

where 

X^'\q) = -^Y.G{k + q)G{k). (22) 

k 

We solve Eqs. (UHll-dlSl) and (EOl)-© self-consistently. 



invariant under the transformation Ci2X ~^ ~Ci2i, which 
transforms Of" to OJ^'^' . In addition, for J = 0, the 
orbital space also has the full rotational symmetry and 
is equivalent to the spin space, and thus all the above 
response functions are the same except for x'^^'^'^^° {q) ■ 

The response functions in the FLEX approximation 
are given by 



^charge 


- ^iXll: 


aiiq) 


+ Xn;22(9)], 


(29) 


x'-'iq) 


^^Xll: 


aiiq) 


+ Xn;22(9)]: 


(30) 


X^^iq) 


= 4[X^2: 


;12iq) 


+ Xj2;2l(9)]: 


(31) 


X^\q) 


= MXl2: 


;12('?) 


-Xj2;2l(9)]: 


(32) 


X'\q) 


= MXII: 


aiiq) 


-Xn;22(9)]: 


(33) 


X^^^^q) 


= MXl2: 


Mq) 


+ X?2;2l(9)], 


(34) 


X^'^'iq) 


= 4[X?2: 


Mq) 


-X?2;2l(9)]: 


(35) 


X^^'^'iq) 


= MXII: 




-Xn;22(9)], 


(36) 



where we have used trivial relations such as Xii ii('?) = 
X22-22('i')- Within the FLEX approximation, we can show 
the relation x^^ (q) = X^"'' (q) even for J 7^ 0, while it 
does not hold in general. 



B. Response functions 



C. Gap equation 



By using obtained x^{q) ^^nd x^{q)i we can calculate 
response functions. The response function corresponding 
to an operator Of^ is given by 

X^(q,za;™) =^ r dre-^^-^+'^--^TrOHr)0^), 
^ Jo 

where o denotes the origin. An operator Of^ in the 
second-quantized form is given by 



The matrix elements O^^.^,^, are given by 



(24) 



/-J charge 
Ta\T' a' 




(25) 


O'^" , , 

^ ra\r' a' 


= ^TT'Ka'^ 


(26) 


0^" , , 


^ tt' ^oa' 1 


(27) 


^ r<7\r' a' 




(28) 



for charge, spin, orbital, and spin-orbital coupled opera- 
tors, respectively, where is the Pauli matrix for v (= x, 
?/, or z) component. Due to the rotational symmetry in 
the spin space, the relations x" {q) — x'^ (?) — x" (?) 
and x^^^'l?) = X^''''\q) = X^^^'l?) hold. In addi- 
tion, the present model has rotational symmetry in the 
T^-T^ plane, and the relations x^ (*?) = X^ (?) ^-nd 
X^ (q) = X^ {q) also hold. We note that there is 
additional symmetry for J = 0. For J = 0, the model is 



Now, we derive a gap equation for superconductivity. 
First, we categorize the anomalous self-energy by sym- 
metry. The anomalous self-energy for a spin-singlet state 
is given by 



.singlet/, N _ 



1, 



The anomalous self-energy for a spin-triplet state is given 
by 



0';;P'^*(fc) = ^[0.T;.'i(fc) + '/'.i;r'T(fc)]- 



(38) 



Due to the rotational symmetry in the spin space, the 
spin-triplet states with (/)*"?'°*(/c), with (^t-|.t-/|(A:), and 
with 0^|.T-/j^(fc) are degenerate. We can categorize su- 
perconducting states further by orbital symmetry. For 
an orbital-parallel-antisymmetric state (orbital- state) , 
the anomalous self-energy is defined by 



T,T',T" (39) 

^ -\[4iik) - 42ik)]- 

For an orbital-parallel-symmetric state (orbital-r^ state) , 
the anomalous self-energy is defined by 



(40) 



4[0ii(fc) + 4(^)]- 
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For an orbital-antiparallel-symmetric state (orbital-r^ 
state), the anomalous self-energy is defined by 



(41) 



For an orbital-antiparallel-antisymmetric state (orbital- 
r° state), the anomalous self-energy is defined by 



4{k) = -'^ E <r'Sr'r"4''Ak) 

t^t'.t" 

-^[0f2(fc)"4l(fc)]- 



(42) 



An orbital-r"^ state with = x, y, or 2: is a state with 
a d- vector for the orbital space parallel to the v axis, 
while an orbital-r" state is an orbital-singlet state. The 
anomalous self-energy has the following symmetry: 



0rcr;T'(T'( fc) — ^t'ct' ;TO- (fc) • 

Then the following relations are obtained: 

^singlct(_^) = C"*5l<=*(fc), 
^triplct(_^) = -0^,"P'^*(fc), 

0*"P'"*(-fc) = 0o"P'°'(fc), 



(43) 
(44) 

(45) 
(46) 
(47) 



where 1/ = x, y, or z. 

The linearized gap equation for the anomalous self- 
energy is expressed as 



A(r,c,j>)4(fc) 



k' 

T 



= -j^J2vHk-k')\G{k')\^4ik'), 

k' 



(48) 



with A(r, ^, !>) ~ 1, where T denotes a representation 
of tetragonal symmetry G^y which 4ik) obeys, D — x, 
y, z, or 0, and F^{k) is defined by the same way as 
4{k)- Thus, the superconducting transition tempera- 
ture is given by the temperature where an eigenvalue 
A(r,^, i>) of Eq. (|48|) becomes unity. The effective pair- 
ing interactions V~{q) are written as 



where 



^i('Z)=^llll('Z)-l4;2l(9), 
Kf('?)=nlll('Z)+nl;2l('?), 
K^('?)=nl22('Z)+ni;12('?), 



(49) 
(50) 
(51) 
(52) 



(53) 



(54) 



Due to the rotational symmetry of the model in the 
T^-T^ plane, the orbital-r^ state and the orbital-r^ state 
are degenerate. In addition, for J — 0, the orbital 
space has full rotational symmetry, and orbital- r^, -t^, 
and -T^ states are degenerate, that is, they are orbital- 
triplet states. Moreover, in this case, the spin space and 
the orbital space are equivalent, and thus, a spin-singlet 
orbital-triplet state and a spin-triplet orbital-singlet state 
are degenerate. Note also that, by changing phases of 
wavefunctions, we can show that a spin-singlet orbital- 
singlet state and a spin-triplet orbital-triplet state are 
degenerate for J — 0. 

Before presenting calculated results, here we briefly 
discuss possible candidates for pairing symmetry of su- 
perconductivity. For s-wave pairing, a spin-triplet state 
is favorable since a solution 01 (fc) does not have to change 
its sign in the k space for a spin-triplet state due to the 
sign in Eq. ((54|) when fluctuations x^{l) and/or x'^il) 
are large. In single-orbital models, such an s-wave spin- 
triplet state is odd in frequency, which hardly appears 
in ordinary cases. However, in the present model, an s- 
wave spin-triplet state with even frequency dependence 
is allowed for an orbital-r" state [see Eq. (|T7)) ]. For a p- 
wave state, a spin-triplet state is unfavorable, since 0|(fc) 
should change its sign in the k space for a p-wave state. 
Thus, a spin-singlet state is favorable for p-wave pairing, 
which is allowed in the present model for an orbital-r" 
state with even frequency dependence [see Eq. (l46l) ]. For 
a d^2_y2 state, a spin-singlet state is favorable because 
of similar logic for p-wave pairing, while such a state is 
an orbital-T^, -r^, or -r^ state. 



IV. RESULTS 

In this section, we show results for a 64 x 64 lattice. In 
the calculation, we use 2048 Matsubara frequencies. In 
this study, we fix the value of the intra-orbital Coulomb 
interaction U = 6t and vary J (= J'). Then the inter- 
orbital Coulomb interaction is given by U' = U — 2J. 
In the foUowings, we discuss superconducting states with 
even frequency dependence only, since we find that eigen- 
values A(r,^, i>) for odd-frequency states are small. 

Figure [T] shows static susceptibilities x^(q) = 
X'^(q, iujm = 0) for J = 0,t, and 2ta.tT = O.OObt and the 
electron number n = {Ntot)/N — 1 per site. Among the 
susceptibilities, the spin susceptibility x'^ (l) is strongly 
enhanced by increasing J, that is, such magnetic fluctu- 
ations are enhanced by the Hund's rule coupling. On the 
other hand, (q) and (q)j which include orbital 

fluctuations, are suppressed by the Hund's rule coupling. 
The charge susceptibility x'^^^'^^'^ i'i) is enhanced a little 
by the Hund's rule coupling, but it is still small. Thus, 
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FIG. 1: q dependence of the susceptibilities for J = 0, t, and 
2t at T = 0.005t, n = 1, and (7 = 6t. 



among various fluctuations, the spin fluctuations for a 
large J are important in the present model. 

In Fig. El^a), we show n dependence of x'^ (Qmax) for 
J = 0, t, and 2i at T = 0.005i, where qmax is defined 
as the wavevector at which x'^*(q) becomes the largest. 
For comparison, we also show x'^ (Qmax) for the non- 
interacting system. The spin susceptibility is enhanced 
by the Coulomb interaction, and is further increased by 
the Hund's rule coupling as is already shown in Fig. [1] 
for n — I. However, qmax does not depend so much on 
the Coulomb interaction and the Hund's rule coupling as 
shown in Fig. [IJb). This fact indicates that the char- 
acteristic wavevector within the FLEX approximation is 
almost determined by the property of the non-interacting 
system, i.e., the Fermi-surface structure. As shown in 
Fig-IHb), the characteristic wavevector is qmax — (ti'iTt) 
around n — 2, and becomes qmax — (tt, 0) by decreasing 
n to n ~ 0.8. Thus, it may be expected that a supercon- 
ducting state with d.j.2_y2 symmetry appears for a large 
n and a p-wave superconducting state appears for a small 
n. We also expect an occurrence of an s-wave state if the 
enhanced spin-fluctuations in a large part of the q space 
are available. 

Concerning the orbital state, only the orbital- r'' states 
are possible for the s-wave spin-triplet and p-wave spin- 
singlet states as is discussed in the previous section. For 
d2,2_j,2-wave spin-singlet pairing, there are three possible 
orbital states, r^, r^, and t^. By increasing the Hund's 
rule coupling, spin fluctuations become dominant, and 
superconductivity is mainly determined by Xii-ii(9) o^' 
Xii;22(9) [see Eq. For J = t and 2t, from the 

calculated results shown in Fig. \T\ xj^" (q) < x^ (q) 
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FIG. 2: (a) Spin susceptibilities x"^ (q) at q = qmax and (b) 
qmax for J = 0, t, and 2t as functions of n at T = 0.005t and 
U = 6t. The thin solid lines represent those for (7 = 0. 



and X^^'^' {^) < X'^''^' i'i), we obtain 

Xll;ll(q) - Xl2;2l(q) < Xn;ll(q) + X'l2;2l(q), (55) 
Xn;22(q) + Xl2;12(q) < Xliai(q) + Xl2;2l(q): (56) 

respectively. Thus, from Eqs. (|i^ - ([?T|) . we expect that 
an orbital-T^ state is the most favorable state among 
the (i^2_j,2-wave spin-singlet states. However, the dif- 
ference among V^™^^'^^{q) [v = x, y, and z) would be 
small, since only Xii;ii(q) and Xii;22(q) ^-^e large and 

X'"'(q) > X^''"'(q) > means x!i;ii(q) - x!i;22(q)- 

Figure [3] shows q dependence of the effective pairing 
interactions Vp (q) = (q, iujm = 0) at zero frequency 
for J = 0,t, and 2t ai T = O.OOSt and n = 1. The mag- 
nitude of the effective interactions are increased by the 
Hund's rule couphng. As is discussed above, VjJ'"^'''*(q) 
is shghtly larger than y|'"gi<'*(q) and l/"'"siet (q) ^ Yot 
J — 0, Vq"^^^'^^ (q) is repulsive for s-wave pairing, but 
becomes attractive for J = 2t. For J = 0, all the suscep- 
tibilities are the same except for the charge one, and we 
obtain 



Xll;22(q) -Xl2:12(q) 



-Xl2;12(q) 

-x"'(q)/4<0. 



(57) 



From the inequalities x^^^^^^i'i) < X^ (q) *md < 
X^"(q) = 4xi2a2(q). we obtain Xn;22(q) < < 
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FIG. 3: q dependence of the effective pairing interactions 
V/(q) for J = 0, t, and 2t at T = 0.005t, n = 1, and (7 = 6t. 



X^2;i2(q)> that is, 



Xll;22(q)-X'l2;12(q) <0, 



(58) 



for J ~ 0. Thus, Vjj*"'''''* (q) is repulsive for J = [see 
Eqs. dnH) and dM])]. For J = 2t, Xii;22(q) is large, and 
the spin-triplet orbital-r*^ pairing interaction becomes at- 
tractive. 

Among all the possible superconducting states, we 
find that only three states, the s-wave (Aig symme- 
try) spin-triplet orbital-r" state, the p-wave (E„ sym- 
metry) spin-singlet orbital- r° state, and the dj.2_j^2-wave 
(Big symmetry) spin-singlet orbital-r^ state, have the 
largest eigenvalues of the gap equation Eq. with 
A(r,^, !>) > 0.5 for some parameter sets we have used. 
Figures ma)-[Hc) show n dependence of eigenvalues for 
these states and x"^* (Qmax) at T = 0.005i for J = 0, t, 
and 2t, respectively. The eigenvalues are enhanced by the 
Hund's rule coupling as the spin susceptibility. In par- 
ticular, for J = 2t, A(Aig, triplet, 0) and A(E„, singlet, 0) 
become larger than unity, that is, superconducting tran- 
sitions occur for these symmetry with transition temper- 
atures higher than 0.005t. The p-wave superconductiv- 
ity appears around n = 0.9, where the characteristic 
wavevector is Qmax — (tt, 0) [see Fig. ^b)]. For the s- 
wave state, n dependence of A(Aig, triplet, 0) is rather 
moderate, since the s-wave pairing utilizes fluctuations 
in a larger part of the q space and it is relatively in- 
sensitive to qmax- By increasing n, the magnitude of 
fluctuations are enhanced, and A(Aig, triplet, 0) becomes 
large. As shown in Fig. [Jc), A(Aig, triplet, 0) is still 
large for small n, and the s-wave state may realize in a 
wider n region if we lower a temperature. The eigenvalue 
A(Big, singlet, y) for the d^^-y^-wave state also becomes 
large by increasing n, however, the spin susceptibility en- 
hances more rapidly. Here, we define a magnetic transi- 
tion temperature where x"^ (Qmax) becomes 50 /t. Then, 
the (i2;2„j^2-wave state does not appear within the calcu- 
lated parameter region, but it may appear by lowering 
temperature and/or by adjusting parameter J. 

Figure inila) shows the highest transition temperatures 
among all the possible superconducting and antiferro- 
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FIG. 4: Eigenvalues A (Aig, triplet, 0), A(Big, singlet, y), and 
A(Eu, singlet, 0), and the spin susceptibilities (qmax) as 
functions of n for (a) J — 0, (b) J = t, and (c) J = 2t at 
T = O.OOSt and U = 6t. 



magnetic ones as functions of n for J = 0, i, and 2t. 
Figure [5Kb) shows the highest transition temperatures 
for J — 2t around n — 1. The antiferromagnetic phase 
extends by increasing the Hund's rule coupling. For 
J = and t, we cannot find any superconducting phase 
within T > 0.005t. For J = 2t around n = 0.9 where 
Qmax — (tt, 0), thc p-wavc spin-singlct orbital-r'^ state ap- 
pears. The superconducting transition temperature for 
the p-wave state probably decreases rapidly at n < 0.8 as 
is expected from n dependence of A(E„, singlet, 0) shown 
in Fig. HJc), while we cannot obtain reliable estimation 
of transition temperatures with low values at present 
because of computational limitations, in particular, a 
lattice-size limitation. For J = 2t around n ~ 1.1 where 
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TABLE I: Equivalent pairing symmetry among s-orbital, p- 
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FIG. 5: (a) Transition temperatures for U = 6t. Solid cir- 
cles, triangles, and squares denote antiferromagnetic transi- 
tion temperatures for J — 0, t, and 2t, respectively. Open 
circles and triangles denote superconducting transition tem- 
peratures for the s-wave (Aig) spin-triplet orbital-r" state and 
for the p- wave (E„) spin-singlet orbital-r" state, respectively, 
for J = 2t. (b) Transition temperatures around n = 1 for 
J = 2t. 



the antiferromagnetic transition temperature tends to 
zero, the s-wave spin-triplet orbital-r*^ state appears. 
As is expected from n dependence of A(Aig, triplet, 0) 
shown in Fig. [Jfc), the superconducting transition tem- 
perature for the s-wave state depends slightly on n. The 
s-wave state may appear also at n < 0.8 with low tran- 
sition temperatures as is expected from the behavior of 
A(Aig, triplet, 0) and A (E„, singlet, 0). 

We have determined the highest transition tempera- 
tures among all the possible superconducting and antifer- 
romagentic ones. To discuss competition and coexistence 
of possible states, we have to compare thermodynamic 
potentials in ordered states, and it is one of future tasks. 



V. PAIRING SYMMETRY IN OTHER 
EQUIVALENT MODELS 

In this section, we discuss pairing symmetry in other 
models which are described by Eq. ([1]). So far, we have 
assumed that the orbitals have s-orbital symmetry (s- 
orbital model), but we can also consider models with 
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different orbital symmetry with the same Hamiltonian 
Eq. (HD. 

The model for doubly degenerate Px and py orbitals 
with (ppa) — (ppTr) — t (p-orbital model) is given by 
Eq. ([1]), but the symmetry of the wavefunctions of or- 
bitals is different from that in the s-orbital model. For 
example, by the symmetry operation x — > — x, the sign of 
the wavefunction for the px orbital changes, and then, the 
anomalous self-energy for orbital-antiparallel (orbital-r^ 
and -T^) states transforms to that transformed in the s- 
orbital model multiplied by (—1) under this symmetry 
operation. Thus, the symmetry of the anomalous self- 
energy in the p-orbital model is different from that in the 
s-orbital model. Note that the model with dzx and dyz 
orbitals {t2g orbitals) with (ddir) = (ddS) — t is equiv- 
alent to the p-orbital model including the symmetry of 
orbitals within the square lattice. The model with de- 
generate eg orbitals with (dda) = (ddS) = t (e^-orbital 
model) is also the same model as the s-orbital model ex- 
cept for orbital symmetry. 

In Table [H we show equivalent pairing symmetry 
among s-orbital, p-orbital, and e^-orbital models, e.g., 
the orbital- r° state with A2g symmetry (g-wave) in the 
p-orbital model has the same transition temperature as 
the orbital-r'^ state with Aig symmetry in the s-orbital 
model. However, the wavenumber dependence of the 
anomalous self-energy is the same among the equiva- 
lent states. From the calculated results for the s-orbital 
model, we find that the g-wa.ve spin-triplet orbital-r'^ 
state appears in the p-orbital model, and the dx^-y^- 
wave spin-triplet orbital- r° state appears in the eg-orbital 
model. The p-wave spin-singlet orbital-r" state appears 
both in the p-orbital and the Cg-orbital models. 

Here we discuss another model for s orbitals with 
two Fermi-surfaces. Under the transformation Ci2cr 
e^^'^^Ci2a with Q = (tt, tt), or equivalently Ck2cr ^ 
Ck-i-Q2cr, the dispersion in Eq. ^ changes as ek2 ^ 
ek+Q2 — — ek2, while the onsite terms in Eq. ([T]) are in- 
variant. Thus, the s-orbital model changes to the model 
described by Eq. ([T]) with the dispersion eti = ^ek2 = 
2t(coskx + cos ky) by this transformation. 

In such a two- Fermi-surface system, we expect super- 
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FIG. 6; Fermi surfaces (solid curves) for the dispersion tki = 
~^k2 = 2t{coskx + cosky) for n = 1. A pair carrying the 
total momentum Q is composed of electrons denoted by closed 
circles on the different Fermi surfaces. 



unchanged. 

Thus, from the calculated results for the s-orbital 
model, we obtain s-wave spin-triplet and p-wave spin- 
singlet states with the finite total momentum Q, even 
without a magnetic field, in the model with two Fermi- 
surfaces which locate around k = (0, 0) and k = Q. At 
temperatures lower than the transition temperatures, a 
superconducting state may coexist with an antiferromag- 
netic (spin density wave) state. A superconducting state 
with the total momentum Q may be stabilized further 
by coexisting with a spin density wave with Q)^i2i and 
it is an interesting future problem. 

Note also that a model for and Py orbitals with 
— {ppa) =^ (ppT^) — t, in which Fermi surfaces locate 
around k = Qi = (tt, 0) and k = Q2 = (0,7r), is equiv- 
alent to the p-orbital model with (ppcr) = (ppTr). We 
can show this equivalence by applying the transformation 
Ciia — > 6**^1 ■'■'CiiCT and Ci2a — > e''^2 '"*Ci2cr. Thus, we also 
obtain superconducting states with the total momentum 
Q = Qi + Q2 in the model with —(ppcr) = [ppi:) — t. 
This conclusion also applies to a model for t^x and ty^ 
orbitals with {ddir) — —{dd6) — t. 



conducting states with the finite total momentum Q like 
the FFLO state by the following reason. In ordinary 
cases, such a superconducting state with a finite total 
momentum q is hard to be realized, since the electron 
on the Fermi surface with the wavevector k has to pair 
with the electron with — k-|-q, but it is off the Fermi sur- 
face unless k satisfies conditions depending on q and the 
structure of the Fermi surface. On the other hand in the 
model with two Fermi-surfaces, when the electron with k 
is on a Fermi surface, the electron with — k-|- Q is always 
on the other Fermi surface as shown in Fig. [5) Thus, 
superconducting states with the total momentum Q is 
naturally expected in the model with two Fermi-surfaces, 
where Q connects centers of these Fermi surfaces. 

Indeed, under the transformation Ck2(T — > Ck+Q2o-, an 
anomalous Green's function for a zero-total-momentum 
state changes as 

-Fl(T;2a' (k, t) = -(T^CklCT(T)c_k2(T') ^^^^ 
FiCT;2a'(k,T; Q) = -(T^CklCT(T)c_k+Q2cr'), 

where Fict;2ct' (k, t; Q) denotes the anomalous Green's 
function for a pairing state with the total momentum 
Q. On the other hand, Fi„.i„i{k,T) and F2(t;2ct' (k, t) 
change to Fi^-i^i (k, t) and F2a-2(T' (k-l-Q, r), respectively, 
i.e., these anomalous Green's functions are correspond- 
ing to pairing states with zero total momentum. Thus, 
by this transformation, the orbital-antiparallel supercon- 
ducting states (orbital-r^ and -t^ states) change to those 
with the finite total momentum Q, while the orbital- 
parallel states (orbital- and -r*' states) remain those 
with zero total momentum. Concerning the susceptibili- 
ties, x^(<l) transforms to X^(q+Q) for A — r^, r^, r^a^, 
and T^a'^ , while the other susceptibilities are unchanged. 
Then, the antiferromagnetic transition temperatures are 



VI. SUMMARY 

We have studied a two-orbital Hubbard model on a 
square lattice. In this model, we have considered orbitals 
with s-orbital symmetry. In such a multi-orbital sys- 
tem, even-parity spin-triplet states and odd-parity spin- 
singlet states are allowed even within states with even 
frequency dependence. Indeed, we have found that the 
s-wave spin-triplet orbital-antisymmetric state and the 
p-wave spin-singlet orbital-antisymmetric state naturally 
appear within a fluctuation exchange approximation. 

Tendencies toward a s-wave spin-triplet state have 
been found also in similar two-orbital models on infi- 
nite dimensional lattices by using a dynamical mean-field 
theoryii^ii& This fact indicates that tendencies toward s- 
wave spin-triplet states are common to two-orbital mod- 
els irrespective of dimensionality. 

We have also discussed pairing symmetry in equivalent 
models which can be described by the two-orbital Hub- 
bard model. The equivalent models with orbital symme- 
try other than s-orbital one require special conditions, 
e.g., ippa) — [ppir) for the p-orbital model, and may 
be unrealistic. However, we believe that these models 
provide an interesting view point to discuss and deter- 
mine pairing symmetry in multi-orbital systems. Note 
also that, in a realistic model such as a model for p or- 
bitals with (ppa) 7^ (pptt), orbital states of a pair mix in 
general, and then, k dependence of the anomalous self- 
energy is not that of an irreducible representation of a 
system. 

Finally, we have discussed models with two Fermi- 
surfaces. In such multi-Fermi-surface systems, we nat- 
urally expect a superconducting state with a finite total 
momentum like the FFLO state. Indeed, we have shown 
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that the superconducting states found in the s-orbital 
model are corresponding to pairing states with a finite to- 
tal momentum, even without a magnetic field, in a model 
with two Fermi-surfaces. Although such a system with 
multi-Fermi-surfaces with the same structure is unreal- 
istic, we expect that such exotic states will be realized 
also in realistic models with multi-Fermi-surfaces whose 
centers locate around different k points with similar, but 
not exactly the same, structures. 
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